Interactions and Disorder in Quantum Dots: Instabilities and Phase Transitions 
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Using a fermionic renormalization group approach we analyse a model where the electrons diffus- 
ing on a quantum dot interact via Fermi-liquid interactions. Describing the single-particle states 
by Random Matrix Theory, we find that interactions can induce phase transitions (or crossovers 
for finite systems) to regimes where fluctuations and collective effects dominate at low energies. 
Implications for experiments and numerical work on quantum dots are discussed. 
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A variety of finite many-body physical systems are 
so complex that one can only approach them statisti- 
cally. Random Matrix Theory (RMT) §] has emerged 
as a unifying language to describe nucleii, atoms, and 
molecules [D at intermediate energies. In these systems, 
the low-energy physics is well described in terms of single- 
particle levels forming shells, but many-body states at 
high enough excitation energy are randomized by inter- 
actions 1^, and well described by RMT. More recently, 
RMT has been applied to quantum chaos and meso- 
scopic physics, including quantum dots [§J^. Here, in 
contrast to earlier applications, RMT is a valid descrip- 
tion of single-particle states, with the external disorder 
potential now inducing the randomization between these 
states. What happens when disorder and interactions are 
both present? This question is relevant to the tunnelling 
of electrons through realistic quantum dots, which has 
been the focus of much recent experimental ||^,^ and the- 
oretical [^J^J^-^ interest. In this Letter we show, using 
renormalization group (RG) methods, that the introduc- 
tion of interactions into quantum dots can produce phase 
transitions (actually sharp crossovers for finite system 
size) in the limit of weak disorder, leading to behavior 
qualitatively different from the noninteracting case. 

Let us focus on a two-dimensional quantum dot (QD) 
of linear size L. The average single-particle level spacing 
is A, the mean free path is Imf, and the diffusion constant 
is D — vplmf/'2, where vp — kp/m* is the Fermi veloc- 
ity, kp \s the Fermi momentum, m* is the effective mass 
of the electron, and fi has been set to 1. The diffusion 
constant is related to another important energy scale, 
the Thouless energy Et = Djl?, which is the inverse 
of the diffusion time through the QD. When the QD is 
strongly coupled to leads the incoming electron's energy 
becomes uncertain upto Et, and therefore the electron 
"samples" g — Et/ ^ single-particle states Q, leading to 
a timncUing conductance of ge? jh. 

We will be interested in QDs with weak disorder 
(ff ^ 1), very weakly coupled to the leads. In this case 
|lO| within any energy band of width Ep (the Thouless 
band) the energy levels are coupled by correlations de- 
scribed by RMT. These correlations make the energies 
rigid so that the probability of finding two close levels 
vanishes as \e\ — €2]^ ■ Depending on the symmetries of 



the problem, there are ten universal random matrix en- 
sembles . We focus on the orthogonal ensemble with 
spin degeneracy for which /3 = 1. Experiments can 
indirectly measure the energy to add an electron to the 
dot (called A2). For the case with spin degeneracy the 
distribution of A2 is the sum of the distribution of sin- 
gle particle spacings (given approximately by the famous 
Wigner surmise) and a delta function at zero energy, to 
account for an opposite spin electron going into a singly 
occupied state, leading to the prediction 



Prmt{^2) = 7;[H^2) + -A2exp(-7rA2)] 



(1) 



This distribution is biniodal and has a peak at zero 
energy. What is observed in experiments has no 
bimodality, and looks more like a symmetric gaussian 
around some nonzero average. These features are also 
observed in numerics which makes it clear that 

the change of shape and shift of the distribution is the 
result of electron-electron interactions. 

The simplest extension of the noninteracting model 
is the constant exchange and interaction (GEI) model, 
where the direct interaction is modelled by a capacitive 
term Q^/2C, and the exchange by a —JS^ term in the 
Hamiltonian. The CEI model (or "universal Hamilto- 
nian" |^,|l^,|l3| Hu) is a very natural starting point in the 
g 00 limit. Using Hu does improve comparison with 
experiments , but the bimodality remains at very low 
temperatures. This does not agree with numerics which 
are carried out at T = and do not show bimodality. 

Motivated by these questions, we have carried out a 
systematic study of how interactions affect low-energy 
properties in a disordered QD. We use a fermionic RG 
method developed by one of us based on the for- 
malism of Shankar , in which one integrates out high- 
energy states and looks at the flow of couplings in the 
low-energy theory. A divergent flow signals a transition 
to a phase whose low-energy physics is not perturba- 
tively connected to the noninteracting limit. For a fi- 
nite QD the phase transition will be replaced by a sharp 
crossover. The advantage of RG over perturbative meth- 
ods is that the presence of regimes of qualitatively differ- 
ent low-energy behavior is made manifest. 

We will assume that the interactions are strong com- 
pared to the disorder. The renormalization of the in- 
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teractions proceeds as in the clean limit leading 
to a Fermi liquid (FL) theory near the Fermi surface. 
For spinless fermions the FL is charaterized by a single- 
particle energy eo(/c) « vp\k — kpl, and an interaction 
function /(fc, k') which parameterizes how the energy of 
a quasiparticle (QP) of momentum k depends on the oc- 
cupation of other QPs 



The only matrix elements with nonzero average are < 



e{k)^eoik) + J2fik,k')dn{k') 



(2) 



with 5n{k) denoting the change of occupation of state k 
with respect to the filled Fermi sea, which is the ground 
state. The low-energy effective FL Hamiltonian is 

Hfl = eo(fc)"(fc) + IY1 /(^' ■ Hk)n{k') ■ (3) 

k kk' 

where n{k) = c^{k)c{k), and the creation and annihila- 
tion operators obey canonical anticommutation relations. 
We will focus on the Thouless band centered on the Fermi 
surface. In particular, there are g momemtum states in 
this band which are completely and randomly mixed ||ic| ] 
to produce the eigenstates of the singlc-particlc Hamil- 
tonian (including disorder) 4>a(k), with energies Ca- In 
terms of these states the disordered FL Hamiltonian is 



HdFL = ^ eaV'lV'a + ^ X! VafljStpi^pltpjtpS (4) 



We de-dimensionalize the FL interaction by f{k,k') — 
ETu{6,d')/g. Note that u depends only on the angular 
position. One resolves the interaction into Fourier com- 
ponents 

oo 

u{0, 9') = u{e - 9') =uo+Y^ u,n cos{m{9 - 6')) (5) 

m— 1 

We will be interested in the case with spin rotation in- 
variance, where there are two FL functions, /s for the 
singlet channel, and ft for the triplet channel. The cor- 
responding (anti)symmetrized matrix elements are 

VaA = T E u,AO,0'){<f>l{k)rp{k') ± cj>l{k')<f>}{k)) 

k,k' 

{^s{k)Mk')±^s{k')cf>-,(k)) (6) 

where the -I- sign is to be taken with the singlet and the 
— with the triplet. Parameterizing the FL interaction as 
@ f{k,k',a,a') = $(fc,fc') + a ■ a'Z{k,k'), where the 
spin of the electron is S* = ^a, the singlet and triplet 
interactions can be expressed |l8j as — ^ — 2>Z and 
-uj = $ + Z. 

In RMT, wavefunctions at different energies are uncor- 
related (except for orthogonality). Denoting an ensemble 
average by <> we have, for the orthogonal ensemble 



< (j>l{k)Mk') >=< M~k)(t>p{k') >-- 



(7) 



UsoA, < V^%^ >= utoA, and < > 



(^) 



Us{'k)1^/ g. The first two terms when rearranged produce 
the CEI model while the last term is a Cooper inter- 
action, which however, is suppressed for our Fermi liquid 
by Since g is assumed large we ignore it henceforth 
§. 

One can similarly calculate the variances of the matrix 
elements. Here one has to pay attention to the subtle 
correlations between different eigenf unctions induced by 
orthogonality . One finds that Uso and ujo produce no 
fluctuations. The general expression is 
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Kurland et al argued [|lj| that since these fluctuations 
vanish in the limit g ^ oo the universal Hamiltonian is 
exact in this limit. On the other hand, because of their 
large number, fluctuation terms can have important ef- 
fects |17 . The best way to see whether the fluctuation 



terms are important for low-energy physics (such as the 
conductance peak spacing experiments) is to carry out 
an RG analysis and see how these terms scale. If in- 
teractions beyond the universal Hamiltonian grow, they 
will dominate the low-energy physics no matter how weak 
they are in the microscopic Hamiltonian. 

We begin with g states in the Thouless band around 
the Fermi surface. We successively integrate out single- 
particle states farthest from the Fermi surface, and at 
some stage in the procedure we have g' states. We will 
define the flow parameter I of the RG by g'{l) = ge^^ 
(with 21 = '^'^ rescale the energy, since 

that would also rescale A, which we want to retain as a 
physical parameter. We want to integrate out the states 
at ±g'/2 to 1-loop order. However, all 1-loop diagrams 
have two internal lines. One of these lines is the state 
at ±5'/2. The other should be summed over all lower 
energy states \g"\ < |5'|/2. This abolishes all reference 
to the states g' /2 and higher in the effective theory. Each 
set of internal lines contributes the following change to 
the interaction matrix elements 



dl 



4 ^"A'l'T ^fjufiS 
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g'G2{n,i^) 



(a^/3) 



(9) 



where V'-''+^^ = V'^"^ + Gi(/^,t/) = (1 - np^ - 

nFu)/{ep. + El.) and G2(^, i^) = {npi, - ?^_F^)/(e^ - e^)- 

In order to make further analytical progress we ensem- 
ble average the internal lines of the diagram (^, v). The 
rationale is that in the course of integrating out states we 
will have to sum over many many states (/u, v). In order 
to carry out the average we use 
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FIG. 1. Phase diagram in the g ^ oo limit. Four different 
phases emerge. 
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At any stage in the RG it is g, and not g' , that appears 
in these averages, since at one-loop there is no wave- 
function renormahzation, and the original wavefunction 
correlations remain unchanged. With this averaging, the 
change in interaction matrix elements can be absorbed 
into a change in the FL parameters. The resulting flow 
equations to zeroth order in 1/g (for n 7^ 0) are 



dl 



dl 



Ce-M(u,„(0)2-|(u,„(/) + Ut„(0) 



Ce-' {Utn{l)f + + Utn{l)f (11) 



where C = log (2) comes from the integration over ener- 
gies less than g' . Each moment n of the Fermi liquid func- 
tion is not mixed with any other moment n', and Uso, Mto 
do not flow. The flow equations have an explicit depen- 
dence on the flow parameter /, but this dependence can 
be removed by defining a new coupling Un{l) = e~'M„(Z), 
which are the FL parameters for a "rescaled" QD with 
a rescaled Thouless energy Elp = e~^ET (corresponding 
to g' = e~^g completely randomized states). Rewriting 
the above equations in terms of $„ = (t2s„ + 3'Sf„)/4 and 
Zn = [utn - '«sn)/4 wc obtain 



d<i>„(l) 

dl 
dZ„{l) 

dl 



-Z„ 



2CZ^ 



(12) 



The decoupling of $„ and Z„ is a consequence of spin 
rotation invariance. From these flow equations we now 
construct the phase diagram that is the central result of 
our paper. Clearly, there are unstable fixed points at 
= — 1/2C and Z„ — — 1/2C which separate weak and 
strong coupling phases. There are four phases (Figure 
1): (i) A phase where both <i>„ and Z„ have small fluc- 
tuations, whose low-energy physics is controlled by the 



universal Hamiltonian. (ii) A phase where Zn diverges 
but <&„ remains small, in which the capacitance will have 
small fluctuations but the exchange interaction will ac- 
quire large fluctuations, (iii) A phase where <&„ diverges 
but Zn remains small, in which there are large capac- 
itance fluctuations, but the exchange energy has small 
fluctuations, (iv) Finally, a phase in which both $„ and 
Zn flow to strong coupling. 

The same flow equations and phase diagram are ob- 
tained in the unitary case if time-reversal symmetry is 
broken by the coupling of an applied magnetic field to 
the orbital sector only, as happens in GaAs dots. 

Let us consider the implications for a real quantum 
dot. Here g « 10 — 100 is typical, which is not partic- 
ularly large. This means that the RG flow should be 
stopped when g' = ge~^ « 1. In the strongly fluctuating 
phases the fluctuations will increase, but the flow will be 
cutoff at the nonuniversal scale I* ~ log {g) and the pu- 
tative phase transition will be replaced by a crossover. 
Nevertheless, we can expect that fluctuations of the ma- 
trix elements play an important role at low energies. 

The peak spacing distribution can be qualitatively un- 
derstood by a two-level example, considered by Levit and 
Orgad 1^. The resulting peak spacing distribution is 
approximately gaussian with a mean position controlled 
by the constant interaction (which is Uso,uto in our lan- 
guage), and a variance given approximately by |^ 



< A2 >^ 4 [< V^22i > - < V1221 



4A2 



(13) 



If bare interaction matrix elements are used in the above 
equation, their contribution is down by 1/g^ compared to 
the single-particle contribution and negligible (assuming 
that \un{l = 0)1 <C 5). In the regime of weak fluctuations 
this remains true, and the width of the distribution is 
controlled mainly by single-particle effects. However, it 
is clear that in the regime of strong fluctuations the vari- 
ance of the interaction matrix elements can dominate the 
single-particle contribution, and thus the width of this 
distribution should increase with interaction strength, as 
has been observed in numerical work . 

The flow equations also have implications for interac- 
tion effects on the width of the quasiparticle levels [ p^ . 
Based on bare matrix elements, the prediction is that 
quasiparticle peaks should be relatively well-defined for 
energies e < A^/g, while they acquire a width larger than 
the single particle spacing at larger energies. If one is well 
into a strongly fluctuating phase, one expects that this 
threshold energy is pushed down and may even become 
zero. In this case the system would not be a Fermi liquid 
and the conductance peaks will have a flnite width even 
at T = (which has been seen experimentally [po|). 

Consider now FL instabilities. In the clean system, the 
FL becomes unstable whenever any of the Z„ < — 1 
|l8| . In the two-level example, the energy to add a parti- 
cle is controlled by Vappa where one of the levels is occu- 
pied and the other empty. In turn this matrix element is 
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the sum of its average (controlled by Uso, uto) and a fluc- 
tuation part (controlled by Ugn, Utn)- If the fluctuations 
overpower the constant interaction, this matrix element 
can become negative, leading to a "bunching" instability 
pl| , p2[ in which more than one particle can be added to 
the QD. The charge/spin density wave instability of the 
clean system when $„, Z„ < — 1 will be modified by dis- 
order and finite-size effects. The characterization of such 
a ground state is an open problem. 

Our approach can be easily generalized to other ensem- 
bles ||l|,|ll| and the results will be reported in forthcoming 
work. It is tempting to speculate about the connection 
with a putative metal-insulator transition in an infinite 
two-dimensional system [^|4j. However, our analysis is 
limited to the zero-dimensional case. Strong fluctuations 
in the matrix elements would have implications for the 
distribution of persistent currents in mesoscopic struc- 
tures . Also relevant is the effect of a weak magnetic 
field, which enhances fiuctuations 1 26 1 , on the phase dia- 
gram. Finally, for realistic QD's, corrections to universal 
RMT behavior (of order l/(?) Q may be important. 

In closing, RMT is a rich field with applications to nu- 
cleii, atoms and molecules, where interactions random- 
ize excited states We have shown that in QD's in- 
teractions and disorder can combine to produce a state 
in which fiuctuations in interaction matrix elements be- 
come large, in which quasiparticles may cease to be well- 
defined at arbitrarily low energies, and the behavior is 
dominated by collective effects. 
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